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During the period when the Sun was intensely active on October-November 
2003, two remarkable solar neutron events were observed by the ground-based 



00 

OS 



neutron monitors. On October 28, 2003, in association with an XI 7. 2 large 

flare, solar neutrons were detected with high statistical significance (6.4<r) by 

the neutron monitor at Tsumeb, Namibia. On November 4, 2003, in association 

with an X28 class flare, relativistic solar neutrons were observed by the neutron 

monitors at Haleakala in Hawaii and Mexico City, and by the solar neutron 

i/-) , telescope at Mauna Kea in Hawaii simultaneously. Clear excesses were observed 
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Haleakala, and 5.2 a for Mexico City. The detector onboard the INTEGRAL 
satellite observed a high flux of hard X-rays and 7-rays at the same time in these 
events. By using the time profiles of the 7-ray lines, we can explain the time 
profile of the neutron monitor. It appears that neutrons were produced at the 
same time as the 7-ray emission. 
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1. Introduction 

Relativistic particles, in particular solar neutrons, give information about ion accel- 
eration in solar flares. Several observations of solar neutrons in solar cycle 23 have been 
reported using the international network of neutron monitors (Usoskin et al. 1997; Lock- 
wood & Debrunner 1999) and solar neutron telescopes (Tsuchiya et al. 2001; Valdes-Galicia 
et al. 2004). Solar flares which produce neutrons frequently occur when activity is near its 
maximum during a solar cycle. In most cases they have also produced X-class solar flares. 
More than one hundred X class flares have been recorded in this solar cycle. 

Intense solar activity occurred from late October to the beginning of November 2003. 
The events that occurred in this period were observed by numerous satellites and detectors 
and have been analyzed by many investigators. During the period when three active regions 
appeared simultaneously on the Sun the soft X-ray flux was very intense and a series of 
eleven X-class solar flares occurred in NOAA regions 10484, 10486 and 10488. At this time, 
solar neutrons were observed on October 28 and November 4, 2003, in association with X17.2 
and X28 class solar flares, respectively. 

An X17.2 solar flare on October 28, 2003, was a remarkable event in this solar cycle. 
Not only was this a large event, but many phenomena were observed in association with 
this flare. It is particularly worth noting the large flux of relativistic particles at the Earth 
(Veselovsky et al. 2004; Panasyuk et al. 2004). Among these particles were solar neutrons 
which were observed by the ground-based neutron monitor before the main Ground Level 
Enhancement (GLE). This solar neutron event has already been reported and discussed by 
Bieber et al. (2005). In this paper, we compare neutron data with 7-ray data observed by 
INTEGRAL satellite and derive the energy spectrum of neutrons using these 7-ray data. 

The second event occurred on November 4, 2003, and solar neutrons were observed by 
NM64 type neutron monitors located at different places, one at Haleakala in Hawaii and the 
other at Mexico City in Mexico. In addition, solar neutrons were also observed by a solar 
neutron telescope located at Mauna Kea in Hawaii. Thus, it is important that a single model 
can explain the data of three detectors to lead an accurate spectrum of solar neutrons from 
the solar neutron event. 

Simultaneous observations of solar neutrons have been made for a few events. In the 
solar event of June 3, 1982, neutrons were observed by a ground-level detector and by 
spacecraft, simultaneously (Chupp et al. 1987). High energy neutrons were detected by the 
IGY type neutron monitor installed at Jungfraujoch, Switzerland, and low energy neutrons 
and high energy 7-rays were observed by the Gamma Ray Spectrometer (GRS) onboard the 
Solar Maximum Mission ( SMM). 
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On May 24, 1990, solar neutrons were observed by the IGY type neutron monitors 
located at Climax and several stations in North America, simultaneously (Shea et al. 1991; 
Debrunner et al. 1997; Muraki & Shibata 1996). On June 4, 1991, solar neutron signals were 
recorded by the neutron monitor and the solar neutron telescope located at Mt. Norikura 
(Muraki et al. 1992; Struminsky et al. 1994). However, in this event, the energy spectrum 
of solar neutrons calculated from the data of these detectors were not self consistent. This 
discrepancy came from the propagation model of solar neutrons in the Earth's atmosphere. 
By using the same propagation model, which is Shibata model (Shibata 1994), nearly the 
same spectrum was obtained (Shibata et al. 1993). In this paper, we report the analysis 
results of the November 4 event using data from the neutron monitors, the solar neutron 
telescope and the spacecraft. Our model can explain data from three detectors consistently. 

2. Solar Neutron Event Associated with an X17.2 Flare on 2003 October 28 



2.1. Observations 
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Fig. 1. — Top panel shows the time profiles of bremsstrahlung 7-rays observed by the IN- 
TEGRAL satellite on October 28, 2003. Second to fourth panel shows the line 7-ray time 
profiles observed by the INTEGRAL satellite on October 28, 2003. The bremsstrahlung 
component has been subtracted. Second panel shows the time profile of the 2.2 MeV neu- 
tron capture 7-rays. Third panel gives the profile of 4.4 MeV 7-ray of C nuclei, and the 
fourth one that of the 6.1 MeV 7-ray of O nuclei. The final panel shows the sum of the data 
in the third and fourth panels. 
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Fig. 2. — Spectra of 7-rays between 1.5 and lOMeV observed by INTEGRAL from 
11:02-03 UT (left) and 11:03-15 UT (right) on October 28, 2003, with background sub- 
tracted. In right figure, clear signals of 2.2, 4.4 and 6.1 MeV 7-rays appear superimposed on 
the bremsstrahlung component. 
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An X17.2 class solar flare occurred at 9:51 UT (this is the time observed at the Earth, 
and same definition here after) on October 28, 2003, located in NOAA active region 10486 
at S16° E08°. From 10:36 to 11:06 UT, an interval which includes the start time of intense 
emission of soft X-rays of the X17.2 flare, the RHESSI satellite was, unfortunately, in the 
South Atlantic Anomaly (SAA). However, intense emission of high energy 7-rays was seen 
in the data after 11:06 UT, indicating that strong particle acceleration occurred during this 
flare. 

On the other hand, large fluxes of hard X-rays and 7-rays were observed by the IN- 
TEGRAL satellite shortly after 11:00 UT. Figure 1 show the bremsstrahlung and line 7-ray 
time profiles from INTEGRAL. In the top panel of Figure 1, two peaks of intense emission 
of bremsstrahlung 7-rays are seen at around 11:03 and 11:05 UT. However, there is only one 
peak (around 11:05 UT) in line 7-ray time profiles as shown in the second to fourth panel in 
Figure 1. This more or less coincides with the second peak in the bremsstrahlung 7-rays. 

Figure 2 shows 7-ray spectra between 11:02 and 11:03 UT, and between 11:03 and 
11:15 UT. From 11:02 to 11:03 UT, when the first peak of bremsstrahlung 7-rays was seen, 
there is no line 7-ray component. 7-ray lines were clearly seen in the 7-ray spectrum after 
11:03 UT, consistent with the line 7-ray time profiles shown in Figure 1. Thus, it appears 
that the time profile of ion and electron acceleration were quite different at this event, and ion 
acceleration either didn't occur or was quite weak during the first peak of the bremsstrahlung 
7-rays. We can assume that the ion acceleration only started after 11:03 UT. 

In Figure 1, note that the 2.2 MeV neutron capture 7-ray line peaks around 11:06 UT 
and has a long decay time. The 4.4 and 6.1 MeV 7-ray lines of de-excited C and O ions 
peak around 11:05 UT, giving about a one minute gap between the two peak times. In 
general, neutron capture 7-rays are delayed from 7-ray lines of de-excited ions, since it takes 
time for high energy neutrons to slow down and be captured by protons (Wang & Ramaty 
1974). Thus, it is evident that solar neutrons were produced at this flare, and were probably 
produced at the same time that 4.4 and 6.1 MeV 7-ray lines were emitted. Hereafter, we 
assume that solar neutrons were produced around 11:05 UT. 

At 11:05 UT on October 28, 2003, the Sun was located over Africa. Among our in- 
ternational network of solar neutron telescopes, Gornergrat in Switzerland and Aragats in 
Armenia had a possibility of observing solar neutrons. On the other hand, Tsumeb observa- 
tory (17.6°E, 19.1°S, 1240 m a.s.l.) was located just under the Sun at this time. The altitude 
of Tsumeb Observatory is a little bit low, however the air mass for the line of sight to the 
Sun was thinner than any other observatory because the zenith angle of the Sun was 9.5°. 
Solar neutrons were clearly observed by the Tsumeb neutron monitor (Bieber et al. 2005). 
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Fig. 3. — Top panel shows the five minute counting rate observed by the Tsumeb neutron 
monitor on October 28, 2003. The solid smooth line is the averaged background and the 
dashed lines are ±lcx from the background. The bottom panel gives the one minute counting 
rate of the Tsumeb neutron monitor (black line) and the time profile of the Lomnicky Stit 
neutron monitor (gray line). The solar neutron event in the Tsumeb data started well before 
the GLE event seen at Lomnicky Stit. 
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The five minute counting rate of the Tsumeb neutron monitor is shown in top panel 
of Figure 3. Clear excesses are seen between 11:05 and 11:15 UT and between 11:20 and 
11:25 UT. The statistical significance of these excesses are 4.8 a for 11:05— 11:10 UT, 4.2 a 
for 11:10 — 11:15 UT and 3.4 a for 11:20— 11:25 UT. The total significance for the ten minutes 
between 11:05 UT and 11:15 UT is 6.4(7. 

At the same time, high energy protons were produced in association with this flare, and 
large ground level enhancements (GLEs) occurred around the world. We can exclude the 
possibility that excesses observed at Tsumeb came from energetic ions, by considering the 
time profile of the neutron monitor at Lomnicky Stit (20.2°E, 49.2°N, 2634 m a.s.l.) together 
with that of the Tsumeb's neutron monitor (bottom panel of Figure 3). The start time of 
the first excess of the Tsumeb neutron monitor is about 10 minutes earlier than the event at 
Lomnicky Stit, while the second excess at Tsumeb is consistent with this event at Lomnicky 
Stit. Thus, it appears that the second excess at Tsumeb came from energetic ions and the 
first excess was solar neutrons. 



2.2. Analysis result 

Using observational data presented in Section 2.1, we can calculate the energy spectrum 
of the solar neutrons even though neutron monitors cannot measure the energy of neutrons. 
By using the time of flight (TOF) method and assuming the emission time of solar neutrons, 
a spectrum can be derived. We assume that the neutrons were produced at 11:05 UT, when 
line 7-ray emission peaked and that the energy of the neutrons responsible for the excesses 
recorded by the neutron monitor is greater than 100 MeV. 

From the time profile of the neutrons, we calculate the energy spectrum of solar neutrons 
at the solar surface by the following formula: 



ex S x AE n x P w 

where AN is the number of excess counts contributed by solar neutrons, and e is the detection 
efficiency of the neutron monitor. Here, e includes the attenuation of solar neutrons through 
the Earth's atmosphere. S is the area of the neutron monitor, AE n is the energy range 
corresponding to one time bin and P is the survival probability of solar neutrons traveling 
from the Sun to the Earth. 

To obtain the e of Equation (1), we calculated the attenuation of solar neutrons by the 
Earth's atmosphere using the Shibata model (Shibata 1994). Solar neutrons with energy 
less than 100 MeV are strongly attenuated by the Earth's atmosphere, so the detection 
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of neutrons by the Tsumeb monitor directly implies that the spectrum extended beyond 
100 MeV. For the detection efficiency of the neutron monitor, we used calculated result by 
Clem & Dorman (2000). 
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Fig. 4. — The energy spectrum of neutrons at the solar surface on October 28, 2003. 
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Using these observational and simulation results, we calculated the energy spectrum of 
neutrons at the solar surface using the method used for the solar neutron event observed on 
November 24, 2000 (Watanabe et al. 2003). Figure 4 shows the result from Equation (1). 
By fitting these data points with a power law of the form C x (E n /100[MeV]) a , the energy 
spectrum of solar neutrons was obtained. The energy spectrum is fitted by a power law as; 



For this fit, xVdof = 7.10/8 = 0.89, and the \ 2 probability is 53 %. The fitting region is 
above 100 MeV. This power index is a typical value for solar neutron events observed thus 
far. The total energy flux of > 100 MeV neutrons emitted by the Sun was estimated to be 
3.1 x 10 25 erg/sr. 



In the analysis method described above, the energy spectrum is calculated by divid- 
ing the response into several bins, each characterized by a mean energy. For the survival 
probability of solar neutrons, as well as the attenuation of neutrons and detection efficiency 
of the neutron monitor, the values at these discrete energies are used. In order to calcu- 
late the energy spectrum of the solar neutrons in detail, we include an assumption about 
the time profiles of solar neutrons but still assume a power law spectral index at the solar 
surface. Using this method, we can investigate whether the neutrons are produced continu- 
ously. To clarify the consistency with the conventional method, we begin by assuming that 
the neutrons are produced impulsively. 





2.2.1. Simulation by Impulsive Model 
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Fig. 5. — Top panel shows the reduced ^-distribution of the fit of the simulated counting 
rate to the observed excess of the Tsumeb neutron monitor. Two minute counting rate is 
used in this calculation. The X-axis represents the power index of the simulated time profiles, 
and Y-axis corresponds to the ratio of simulated counting rate (N) to the observed one (N Q ). 
In this fitting, data obtained during 11:05— 11:15 UT are used. Bottom panel shows the two 
minute counting rate (solid line) observed by the Tsumeb neutron monitor on October 28, 
2003, together with simulated time profile (points) for which solar neutrons are assumed to 
be produced instantaneously at 11:05 UT, when power index is —3.5. 
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In this simulation, the power index of the neutron spectrum at the Sun is changed from 
— 1.5 to —7.0 with a step of 0.1, while the energy range of the incident neutrons confined to 
50 — 1500 MeV. The time profile of neutrons detected by the neutron monitor is calculated 
using the neutron attenuation in the Earth's atmosphere given by the Shibata model (Shibata 
1994) and the detection efficiency of the neutron monitor as calculated by Clem & Dorman 
(2000). The decay of neutrons between the Sun and the Earth is also taken into account. 
The result of this simulation can then be compared with the observational data, normalizing 
the simulated counting rate (N) to the observed excess counting rate (N ). 

The top panel of Figure 5 shows the reduced ^-distribution of the fit of the simulated 
counting rate to the observed excess of the Tsumeb neutron monitor obtained from the 
following formula: 



In this fitting, data obtained from 11:05 to 11:15 UT are used, x 2 nas its smallest value 
when the spectral index is around —3.5. When the spectral index is —3.5, the simulated 
result reproduces the observed result as shown in the bottom panel of Figure 5, where 
X 2 /dof = 6.07/4 = 1.52, which yields the minimum x 2 f° r the simulated time profile. From 
this fitting, the energy spectrum is determined as follows: 



This is comparable to the result obtained using the simpler method shown in Equation 
(2), but the total energy flux of solar neutrons with energy range between 50 — 1500 MeV is 
(9.8+J; 2 ,) x 10 25 erg/sr, about a factor of three higher than the other estimate. This is because 
of the lower cutoff energy of the neutron spectra. 




(3) 




+0.4 




2.2.2. Simulation by Neutron Production with ^-ray Time Profile 
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Fig. 6. — The observed and simulated time profiles of the Tsumeb neutron monitor. Solid 
line is the observed two minute counting rate and points indicate the best fit simulated time 
profile for solar neutrons assumed to be produced with the same time profile as 7-ray lines. 
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We next simulated the neutron time profiles detected at Tsumeb by assuming that 
neutrons were produced with a time spread, since extended production of line 7-ray was 
observed by INTEGRAL. For this calculation, we used the summed time profile of 4.4 and 
6.1 MeV 7-ray as shown in the bottom panel of Figure 1, as the production model of solar 
neutrons. These are the 7-ray lines of Carbon and Oxygen, which indicate the time profile 
of ion acceleration. We used the data observed from 11:02:45 to 11:10:00 UT. The spectral 
index of neutrons at the Sun is varied from —1.1 to —7.0 in steps of 0.1. The energy range 
of the neutrons is again taken to be 50 — 1500 MeV. 

The x 2 f° r the fit were calculated by using data obtained from 11:05 to 11:15 UT. \ 2 
has its smallest value for the spectral index —2.9. The simulated result reproduces the 
observed result most closely when the spectral index is —2.9 as shown in Figure 6. When 
the spectral index is —2.9, xVdof = 2.76/4 = 0.69, which provides the minimum x 2 fc> r the 
simulated time profiles. From this fit, the spectral index is found to be —2.9 ± 0.3. The 
best fit spectral index is harder than the index derived on the assumption that the neutrons 
were produced impulsively, but the total energy flux of the neutrons is now estimated to be 
(6.2^) x 10 25 erg/sr, not very much different from the result for impulsive production. 



3. Simultaneous Observations of Solar Neutrons on 2003 November 4 



3.1. Observations 
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Fig. 7. — Spectrum of 7-rays between 1.5 and 10 MeV observed by INTEGRAL from 19:40 to 
19:50 UT on November 4, 2003, with background subtracted. A signal produced by 2.2 MeV 
7-rays appear superimposed on the bremsstrahlung component. There is weak evidence for 
4-7 MeV 7-ray lines. 
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Fig. 8. — Time profiles of 7-ray lines observed by the INTEGRAL satellite on November 4, 
2003. The bremsstrahlung component has not been subtracted. Top panel shows the time 
profile of the 2.2 MeV neutron capture 7-ray. Second panel gives the 4.4 MeV 7-ray of C 
nuclei, and the third one contains the 6.1 MeV 7-ray of O nuclei. The final panel shows the 
sum of the data in the second and third panels. Although these time profiles contain line 7- 
ray components which indicate the time profile of ion acceleration, the dominant component 
is bremsstrahlung. 
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On November 4, 2003, an X28 class solar flare occurred at 19:29 UT, located in NOAA 
active region 10486 at S19° W83°. This is the largest solar flare on record. At around 
19:42 UT, intense emission of soft X-rays was detected by the GOES satellite such that the 
detection was saturated. After 19:42 UT, intense emission of hard X-rays and 7-rays was 
observed by the INTEGRAL spacecraft. Unfortunately, at this time, the RHESSI spacecraft 
was on the night side of the Earth. Figure 7 shows the energy spectrum of 7-rays observed 
by INTEGRAL. In this event, although the components of the line emission produced by de- 
excited ions, C(4.4MeV) and 0(6.1 MeV) were not prominent, the 2.2 MeV neutron capture 
line can be clearly seen. Intense bremsstrahlung X-rays and 7-rays were also observed. 
Figure 8 shows the time profiles of 7-rays for different energy bins that contain line 7-ray 
components produced as a result of the ion acceleration. There is a delay of the 2.2 MeV 
neutron capture 7-ray emission from that of the line 7-ray components produced by excited 
ions of C and O. We may assume that ion acceleration occurred at the same time as the 
7-ray lines were emitted, although the main component of these 7-rays is bremsstrahlung. 
And we may assume that solar neutrons were produced at the same time. 

At 19:45 UT, the Sun was located between Hawaii and South America. Although the 
Chacaltaya observatory was the best place to observe solar neutrons in our international 
solar neutron telescope network at this time, no data are available because of a data gap. 
Sierra Negra would also have been a good place to observe solar neutrons, but the Mexico 
solar neutron telescope had not started continuous observation at that time. Thus it was 
necessary to examine data from the Hawaii observatory, which was the third closest of the 
observatories able to observe solar neutrons. 

At 19:45 UT, the zenith angle of the Sun was 49.9° at Mauna Kea, and 50.5° at Haleakala. 
The air mass along the line of sight to the Sun was 947g/cm 2 and 1112 g/cm 2 , respectively. 
The other suitable location was Mexico City, where the zenith angle of the Sun was 40.52° 
and the air mass along the line of sight to the Sun was 1026 g/cm 2 . Attenuation of solar 
neutrons by the Earth's atmosphere above these observatories is calculated using the Shibata 
model (Shibata 1994). The Haleakala and the Mexico City observatories have nearly the 
same attenuation, while Mauna Kea is located at the best place to observe solar neutrons. 
However, simultaneous signals were found in both the Haleakala and the Mexico City neutron 
monitors. 
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Fig. 9. — Five minute averages of the counting rate observed by the Haleakala (top) and 
Mexico City (bottom) neutron monitors on November 4, 2003. The solid smooth line is the 
averaged background and the dashed lines are ±lcr from the background. 
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Solar neutrons were observed by the 18NM64 neutron monitor at Haleakala, Hawaii 
(203. 7°E, 20.7°N, 3030 m a.s.L). The top panel of Figure 9 shows the five minute averages 
of the counting rate observed on November 4, 2003. At this time, the sampling interval 
of the Haleakala neutron monitor was 10 seconds. Clear excesses were seen after 19:45 UT, 
continuing for 15 minutes. The statistical significances of these excesses are 4.5 a for 19:46:20 
- 19:51:20 UT, 5.3 a for 19:51:20 - 19:56:20 UT, and 3.1a for 19:56:20 - 20:01:20 UT. The 
total significance for the fifteen minutes between 19:46:20 UT and 20:01:20 UT is 7.5 a. Note 
that this time interval was just taken to get the maximum significance. 

Solar neutrons were also observed by the 6NM64 neutron monitor at Mexico City 
(260. 8°E, 19.33°N, 2274 m a.s.L), as shown in the bottom panel of Figure 9. At this time, the 
sampling interval of the Mexico City neutron monitor was 5 minutes. Clear excesses were 
seen after 19:45 UT, which was the same time as the excesses observed by the Haleakala neu- 
tron monitor. The statistical significances of these excesses are 2.6 a for 19:45— 19:50 UT, 
3.1 a for 19:50-19:55 UT and 3.3 a for 19:55-20:00 UT. The total significance for the fifteen 
minutes between 19:45 UT and 20:00 UT is 5.2 a. 
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Fig. 10. — The five minute averages of the counting rate of PMT_L, PMT_H, and 
layer l_with_anti channels (see details in the text) of the Hawaii solar neutron telescope on 
November 4, 2003. The solid smooth line is the averaged background and the dashed lines 
are ±lcr from the background. 
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One would expect that Mauna Kea (203. 7°E, 19.8°N, 4200 m a.s.l.) should be a better 
place to observe neutrons in this event than Haleakala and Mexico City. This is the location 
of the Hawaii solar neutron telescope with an area of 8 m 2 , constructed from proportional 
counters and plastic scintillators, but only a minimal excess was seen after 19:45 UT in the 
PMT_L, PMT_H, and layer l_with_anti channels in this telescope as shown in Figure 10. The 
PMT_L and PMT_H are channels of scintillation counter which detect neutrons (recoil pro- 
tons), energy thresholds of which are 12 MeV and 20 MeV, respectively. The layer l_with_anti 
is proportional counter channel which is located under the scintillation counters. This ap- 
parent discrepancy between neutron monitors and the solar neutron telescope is discussed 
in the next section (Section 3.2.1) considering the surrounding environment of the detector. 



3.2. Analysis result 

In order to understand the Mauna Kea result, we first use the other observational data 
to estimate the energy spectrum of the solar neutrons. We begin with the data from the 
Haleakala monitor, because it recorded the largest excess with the best time resolution. 
We determine the neutron energy by using the TOF method, assuming that all the solar 
neutrons were produced at 19:45 UT, the peak of the intense emission of high energy 7- 
rays observed by INTEGRAL as shown in Figure 8. Under this assumption, the energy 
of neutrons observed by the Haleakala neutron monitor between 19:51:20 and 19:56:20 UT 
ranged from 59 — 913 MeV. 

To derive the energy spectrum of neutrons at the solar surface from the observed time 
profile by the neutron monitor, the survival probability of neutrons between the Sun and 
the Earth, the attenuation of solar neutrons passing through the Earth's atmosphere and 
the detection efficiency of the neutron monitor must be taken into account. Attenuation 
is calculated using the Shibata model (Shibata 1994), and we used the detection efficiency 
calculated by Clem & Dorman (2000). 
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Fig. 11. — The energy spectrum of neutrons at the solar surface on November 4, 2003, 
calculated from the data of the Haleakala neutron monitor. 
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Using these observational and simulation results, we calculated the energy spectrum of 
neutrons at the solar surface using the same method as Section 2.1. The result is shown in 
Figure 11. This spectrum was derived from two minute averages of the counting rate, with 
the vertical errors shown only statistical. The energy spectrum is well fitted by a power law 
as: 



The fitting region is chosen as 100 MeV and above, because there the errors from neutron 
attenuation in the Earth's atmosphere are small. For this fit, x 2 /dof = 0.92/3 = 0.31 and the 
X 2 probability is 82 %. This spectral index is typical of solar neutron events observed thus 
far. The total energy flux of neutrons emitted from the Sun in the energy range 59 — 913 MeV 
is estimated to be 3.4 x 10 26 erg/sr. 





3.2.1. Simulation Using the Impulsive Model 
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Fig. 12. — The best fit simulated time profiles (dots) when the spectral index is —3.9 for 
Haleakala (top), and —4.3 for Mexico City (bottom), superposed on the observed counting 
rate. The start time of the simulated time profile is 19:45 UT, corresponding to the peak 
time of 7-ray emission. 
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By using the same method of Section 2.2.1, time profiles of solar neutrons assuming their 
spectral index at the solar surface was simulated if the neutrons were produced impulsively. 
We examine the x 2 °f the fit of the simulated counting rate to the observed excess of the 
Haleakala neutron monitor obtained from Equation (3). In this fitting, data obtained from 
19:45 to 20:06 UT are used. x 2 nas hs smallest value when the spectral index is around 
—3.9. When the spectral index is —3.9, x 2 /dof = 10.3/7 = 1.47, which yields the minimum 
X 2 for the simulated time profile (the top figure of Figure 12). From this fitting, the energy 
spectrum is determined as follows: 



This is consistent to the result obtained using the simpler method shown in Equation (5). 
The total energy flux of solar neutrons with the energy range between 50 — 1500 MeV is 
(6.7t° i) x 10 26 erg/sr, about the same order as calculated value from Equation (5) . 

We have done the same analysis for the Mexico City neutron monitor. In this fitting, 
data obtained during 19:45— 20:05 UT are used, x 2 nas its smallest value when the spectral 
index is around —4.3. When the spectral index is —4.3, x 2 /dof = 5.55/3 = 1.85, which 
yields the minimum x 2 for the simulated time profile (the bottom figure of Figure 12). From 
this fitting, the energy spectrum is determined as: 



Although the spectral index is softer than the result of the Haleakala neutron monitor, it is 
consistent with Equation (5). The total energy flux of solar neutrons with energies between 
50 — 1500 MeV is calculated to be (5.4 ±0.7) x 10 26 erg/sr, about the same order as the result 
obtained from the analysis of Haleakala. 

We then simulated the time profile of neutrons which should be observed from the Hawaii 
solar neutron telescope using the energy spectrum of incident neutrons obtained from the 
data of the Haleakala neutron monitor, namely 1.5 x 10 28 x (En/100 [MeV])~ 3 - 9 [/MeV/sr]. 
We did not attempt to derive a spectrum from the data, because excesses of the solar neutron 
telescope are small. 




1+0.1 
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Hawaii SNT : Main Channels (concrete 20cm : neutron) 
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Fig. 13. — Detection efficiencies of the Hawaii solar neutron telescope for neutrons when the 
detector is surrounded by a 20 cm concrete wall. The black lines indicate the PMT_L (solid) 
and PMT_H (dashed) scintillator channels. The gray lines indicate layer channels of layerl 
(solid), layer2 (dashed), and layer3 (dash-dotted) with anti coincidence of the anti counter. 
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The detection efficiency of the Hawaii solar neutron telescope is calculated using Geant3, 
FLUKA— COLOR model. In this calculation, the Hawaii solar neutron telescope is sur- 
rounded by the 20 cm concrete walls, since it is constructed within the building housing 
of the SUBARU telescope. Figure 13 shows the detection efficiencies of the Hawaii solar 
neutron telescope for neutrons and 7-rays. 
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Fig. 14. — The simulated time profile (dots) of five minute counting rate of layerl _with_anti 
channel of the Hawaii solar neutron telescope on November 4, 2003, superposed on 



the observational data. The energy spectrum of the incident neutrons is 1.5 x 10 x 
(En/100 [MeV])~ 3 ' 9 [/MeV/sr], which was obtained from the data of Haleakala neutron mon- 
itor. The start time of this time profile is 19:45 UT, corresponding to the peak time of 7-ray 
emission. 
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The simulated result for the layerl_with_anti channel of the Hawaii solar neutron tele- 
scope, which recorded largest excess, is shown in Figure 14. For this fitting, x 2 /dof = 
1.24/3 = 0.41, so the simulation result is consistent with the observed excess. Because of the 
high counting rate from the non-hadronic component (7-ray, muon and so on) in the Hawaii 
solar neutron telescope, these excesses are not statistically significant, although they corre- 
spond to the same total flux of solar neutrons observed by the Haleakala neutron monitor. 
It is however possible that solar neutrons actually produced the little hump in the data. 



3.2.2. Simulation by Neutron Production Using the ^-ray Profile 
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Fig. 15. — The simulated time profiles (dots) when the spectral index is —3.9 for Haleakala 
(top) and —4.3 for Mexico City (bottom), superposed on the observed counting rate of the 
Haleakala neutron monitor. Points are the simulated time profile for solar neutrons assuming 
that they were produced with the same time profile as the high energy 7-rays shown in Figure 
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Next we simulated the neutron time profiles detected at Haleakala and Mexico City by 
assuming that neutrons were produced with a time spread. Calculated method is same as 
Section 2.2.2 For this calculation, we used the 7-ray time profile observed by the INTEGRAL 
satellite during 19:42—19:48:00 UT as the production time profile of solar neutrons, as shown 
in the bottom panel in Figure 8. The x 2 of the fit between observed and simulated time 
profiles of the Haleakala and the Mexico City neutron monitor were calculated. In this 
fitting, data obtained from 19:42 to 20:06 UT are used for the Haleakala and from 19:45 to 
20:05 UT for the Mexico City neutron monitor. For the Haleakala data, when the power 
index is —3.9 (top figure of Figure 15), xVdof = 10.57/7 = 1.51, giving the minimum 
value among the simulated time profiles. The spectral index is determined to be —3.9 ± 0.2. 
For the Mexico City data, when the power index is —4.3 (bottom figure of Figure 15), 
X 2 /dof = 4.28/3 = 1.43, giving the minimum value among the simulated time profiles and 
the spectral index determined to be — 4.31q5- The best fit spectral indices are same as 
those derived by assuming that the neutrons were produced impulsively. The total energy 
fluxes of neutrons are estimated to be (7.0 ±0.5) x 10 26 erg/sr from the Haleakala data, and 
( 5 - 7 -as) x 10 26 erg/sr from the Mexico City data. 



4. Discussion and Summary 

Relativistic neutrons were detected in association with the X17.2 solar flare on October 
28, 2003, and the X28 solar flare on November 4, 2003. The detection of October 28 event 
was made by the Tsumeb neutron monitor, and the detection of November 4 event was 
made simultaneously by the neutron monitors at Haleakala and Mexico City, and also by 
the solar neutron telescope at Mauna Kea. During these events, intense emissions of high 
energy 7-rays were observed by the INTEGRAL satellite. 

In order to investigate the production time of solar neutrons, we compared the solar 
neutron data with the 7-ray data obtained from INTEGRAL. In October 28 event, 7-ray lines 
from neutron capture and excited ions of C and O nuclei were clearly observed and were quite 
different from the time profile of bremsstrahlung 7-rays. It appears that the time profile of 
electron acceleration was distinctly different from the time profile of ion acceleration. From 
the time profile of neutron capture 7-rays, it seems that high energy neutrons were produced 
with the same time profile as 7-ray lines of de-excited ions. In November 4 event, the time 
profiles of 7-ray lines which we have assumed to represent the time profile of solar neutron 
production, cannot be obtained independently, since the bremsstrahlung component was 
strong and the line 7-ray components were buried in bremsstrahlung. However, from the 
time profile of the 2.2 MeV neutron capture 7-rays, it appears that the time profile of ion 
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acceleration was approximately same as that of bremsstrahlung emissions. Assuming that 
solar neutrons were produced at the time when these 7-rays were emitted, we could explain 
the observed excesses. 

If we assume that solar neutrons were produced impulsively at 11:05 UT on October 
28, and at 19:45 UT on November 4, when the 7-ray lines peak, we can derive the energy 
spectrum of solar neutrons at the solar surface from the neutron monitors as data by using 
Equation (2) and (5), respectively. For November 4 event, in order to examine whether 
all excesses observed by the Haleakala and the Mexico City neutron monitor and the solar 
neutron telescope at Mauna Kea can be expressed by one energy spectrum consistently, and 
for more detailed analysis, we simulated time profiles of solar neutrons for these detectors, 
and compared with observed time profiles. All of the simulation results are consistent with 
Equation (5) within the range of error as shown in Equation (6) and (7). Thus, we could 
explain all observations with consistent spectrum. 

Although we can fit the data by assuming that solar neutrons are produced impulsively, 
it is more natural to assume that solar neutrons are produced continuously over a finite time. 
We therefore modeled the time profiles of solar neutrons by assuming that the neutrons 
were produced with the same time profile as 7-ray lines from excited ions. For October 28 
event, the spectral indices thus derived are same as those derived by assuming neutrons were 
produced impulsively, —3.9 for Haleakala and —4.3 for Mexico City. In other events, the 
spectral indices derived by assuming neutrons are produced continuously tend to be harder 
than those derived by assuming neutrons are produced impulsively. However, in this event, 
since the 7-rays have a symmetric time profile centering around 19:45 UT, there is little 
difference between the two models. For November 4 event, the result was that the index 
obtained using the line 7-ray time profile is clearly harder (—2.9) than that obtained using 
the impulsive model (—3.5). Therefore, for these two events, the observations were explained 
by assuming that solar neutrons were produced with same time profile as 7-ray lines. 

Although different spectral indices are obtained with the different approach for October 
28 event, these spectral indices are all consistent with the indices calculated from line 7-ray 
observation by the RHESSI satellite (Share et al. 2004) and INTEGRAL satellite (Tatischeff 
et al. 2005). 

The spectrum of accelerated ions can be calculated from the neutron spectrum using 
the spectrum of escaping neutrons produced by the accelerated ions (Hua & Lingenfelter 
1987a,b; Hua et al. 2002). From the neutron spectra shown in Equation (2) and (5), the 
number of protons above 30MeV would be about 10 32 /sr under the assumption that there 
is no turnover of the spectrum. This is a typical value for solar neutron events observed thus 
far. 
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